The Lorentz-Dirac force from QED for linear acceleration 
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We investigate the motion of a wave packet of a charged scalar particle linearly accelerated by a 
static potential in quantum electrodynamics. We calculate the expectation value of the position of 
the charged particle after the acceleration to first order in the fine structure constant in the h — > 
limit. We find that the change in the expectation value of the position (the position shift) due to 
radiation reaction agrees exactly with the result obtained using the Lorentz-Dirac force in classical 
electrodynamics. We also point out that the one-loop correction to the potential may contribute to 
the position change in this limit. 
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A charged particle radiates when it is accelerated. The 
resulting change in its energy and momentum is described 
by the Lorentz-Dirac (or Abraham-Lorentz-Dirac) force 
in classical electrodynamics 0,0, 0- (See, e.g., Ref. 3 
for a modern review.) Thus, if a charge e with mass m 
is accelerated by an external 4-force F^ Kt , then its coor- 
dinates x^{t) at the proper time r obey the following 
equation: 



d 2 x» _ ^ 
m dr 2 ext 



where the Lorentz-Dirac 4-force P^pj is given by 
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We have let c = 1 and defined a c = e 2 /47r. Our metric 
is = diag (+1,-1, -1,-1). 

Although there are many ways to derive Eq. Q in 
classical electrodynamics for a point charge (see, e.g., 
Ref. ||) , a natural question one can ask is whether or not 
this equation can arise in the h — ► limit in QED. It was 
found in Ref. |fj (after an initial claim to the contrary) 
that the position of a linearly accelerated charged particle 
in the Lorentz-Dirac theory is reproduced by the h — ► 
limit of the one-photon emission process in QED in the 
non-relativistic approximation. (See, e.g., Refs. 0,0] for 
other approaches to arrive at the Lorentz-Dirac theory 
from QED.) In this Letter we describe the generalization 
of this work to a fully relativistic charged particle. The 
details will be published elsewhere Q. 

Consider a charged particle with charge e and mass 
m moving in the positive z-direction under a potential 
energy V(z). We assume that V(z) = Vq = const, for 
z < — Z\ and V(z) = for — Z 2 < z for some Z\ and 
Z2, both positive constants. Thus, there is non-zero ac- 
celeration only in the interval {—Z\, —Z<i). The external 
4-force in Eq. JJ) generated by this potential is given 
by F* xt = -V'(z)dz/dT and F c 2 xt = -V'(z) dt/dr with 
F^ xt — F^ Kt = 0. The Lorentz-Dirac 4-force can be given 



byF* 
with 
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A dot indicates the derivative with respect to t. We have 
defined 7 = (1 — i 2 ) -1 / 2 as usual. 

Suppose that this particle would be at z = at time 
t = if the Lorentz-Dirac force was absent (i.e. ife = 0). 
The true position at t — 0, which we denote 5z and call 
the position shift, can readily be found to lowest non- 
trivial order in Fld by treating the Lorentz-Dirac force 
as perturbation. (It was proposed in Ref. 01 that one 
should treat the Lorentz-Dirac force as perturbation.) 
The calculation can be facilitated by using the fact that 
the total energy, m dt/dr +V(z), changes by the amount 
of work done by the Lorentz-Dirac force. Thus, 
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Rearranging and integrating, and then changing the or- 
der of integration, we obtain 

where the final velocity is denoted z\q. We will show that 
this position shift is reproduced to order e 2 in QED in 
the h — > limit. 

The Lagrangian density for the corresponding field- 
theoretic model is 

C = (D^D^ - (m/K) 2 <p^ - -F^ , (6) 



where F M „ = d^Ay - 



dvA^ and 



d^ + i(eA^ + V„)/tl. 



The field <p describes a charged scalar particle with mass 
m and charge e, the field A^ is the electromagnetic field, 
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and the function V p = V(z)S i1 q is the external potential 
which accelerates the charged scalar particle. 

Now, let A'(p) be the creation operator for the charged 
scalar particle with momentum p with p = p z > in 
the positive z region such that the corresponding mode 
function <J> p (i,x) is well approximated by the WKB ap- 
proximation: 
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where pj_ = (Px,p y ), Po = (p 2 + pi+m 2 ) 1/2 , xj. = (a;,y) 
and where k p (z) = {[po ~ V(z)} 2 - m 2 - p 2 L } 1/ ' 2 is the z- 
component of the momentum of the classical particle at 
z with e = 0. If the initial state in the interaction picture 
is j4.t(p)|0), then, to order e 2 , it evolves as follows: 
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where the creation operators at(k) for the photons with 
momenta hk in the Feynman gauge satisfy 

[<V(k), at (k')] = -g^(2ir) 3 2hk 6(k - k') . (9) 

Note that the momentum hk of the photon is of order 
h because the wave number k rather than the momen- 
tum has the classical limit. The .F(p) is the forward- 
scattering amplitude coming from the one-loop diagram, 
which will not be evaluated explicitly. The momentum 
P of the charged particle in the final state is deter- 
mined using conservation of the transverse momentum, 
p^ = P^ + hk± , and energy conservation, po = Pq + hko • 
The expectation value of the z-coordinate of the wave 
packet at t = is given in terms of the forward scatter- 
ing amplitude T{p) and the photon emission amplitude 
Afj,(p, k) as we describe next. 

Let the initial state \i) be given by 
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where the function f(p) is sharply peaked about a mo- 
mentum in the positive z-direction with width of order h. 
The operators A^p) are normalized so that the condition 
= 1 leads to 
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Hence, the function /(p) can be regarded as the one- 
particle wave function in the momentum representation. 



This suggests that the position at t = in the absence of 
radiation is given by 
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where d P =d P — dp- (See Appendix B of Ref. p for a 
derivation.) We assume in this Letter that /(p) is real 
for simplicity. Hence, we have (z)o = as required. 

The final state resulting from the initial state \i) can 
be found from Eq. JSJ as 
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where we have defined F(p 
G"(p,k) 



(1 + &(p))f(p) and 

A^p,k)f(p)y^/p^(dp/dP). The factor 



\J Po/po dp/dP arises due to the change of the integra- 
tion variables from p to P. One can regard the function 
F(p) as the one-particle wave function in the zero-photon 
sector in the p-representation and the function G M (p, k) 
as that in the one-photon sector with a photon with mo- 
mentum hk in the P-representation. This suggests that 
the expectation value of the z-coordinate is 
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We will present a justification of this formula else- 
where Q. Since /(p) is real, we have 
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We have let P — > p in the second term because this term 
will turn out to be of order H°. (Note that (Pq,P) = 
(po,p) at this order because the photon energy hko and 
momentum hk are of order h.) Then, using the assump- 
tion that the function /(p) is sharply peaked, we find for 
small h 



-hd p ReF(p) 
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where p = (0, 0, p) is now the momentum where the func- 
tion / has the peak. The first term is the contribution 
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from the one-loop correction to the potential rather than 
from radiation reaction. Thus, the second term is identi- 
fied as the position shift due to radiation reaction. There- 
fore the position shift in quantum field theory is 

r^*( P ,k) d P A„(p,k). (17) 
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The photon-emission amplitude ^(p, k) can be ob- 
tained to leading order in h by using the WKB wave 
function (J2J) as in Ref. 0. The result turns out to be 
identical with the amplitude for a classical charged par- 
ticle with final momentum p = (0,0, p) which passes 
through the spacetime point (t,x,y,z) — (0,0,0,0). If 
the coordinates of this particle are (r) at proper time 
r, then the classical current generated by this charge is 
j^i.x) = e(dX^/dt)8[* - X(t)]. Thus, we find 



-4"(p,k) = - j d A xe lkot - l ^f{t,x) 



dr- 
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Note that the right-hand side depends on p through the 
world line X^{t). We write X^ — (t, 0,0, z) from now 
on. Let 9 be the angle between the z-direction and the 
vector k and define t; = t — n ■ x. = t — zcos9, where 
h = k/fco- Then Eq. I|18|l can be written 
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where k = feo- We have inserted a function x(0 which 
takes the value one in the region with d 2 X fl /d£? ^ and 
decreases to zero smoothly for large |£| to regularize the 
integral. By using Eq. lj2U|) for *4 M (p, k) and Eq. (|T$|) for 
dpA^ip, k) in Eq. we find 
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where dtt is the integral over the solid angle for k. (The 
p-derivative is taken with £ fixed.) Thus, the position 
shift Szq is independent of the function x(0- 

In order to show that the quantum position shift Szq 
agrees with its classical counterpart, we first note by us- 
ing d£ jdt — 1 — z cos 9 that 
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Next we obtain, by interchanging the order of differenti- 
ation, 
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Then by the chain rule we find 

dpj i \dpj i \dzj i \d P/( 
These equations can be used to find 
d 2 X^ d fdX^ 



(25) 
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By substituting this expression in the quantum position 
shift (|21fl and changing the integration variable from £ 
to t, we find 



Oc_ 

An 



dt- 



dt 



(1 — z cos 9) 3 

where we have changed the integration range from 
(— oo, +oo) to (—oo,0] because z = for (0, +oo), and 
then integrated by parts. 

By substituting dt = d£ + cos 9 dz in the relation dz — 
zdt + (dz/dp)t dp and solving for dz, we obtain 
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dz\ 



(28) 
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By using this expression in Eq. I|27|l and performing the 
51-integration, we obtain 
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where Fld is defined by Eq. (J3J). 

In order to interpret this formula let us consider the 
following equation: 



d z „ , i dt 

m— = [F cxt (z, t) + F LD ] — . 



dr 2 



dr 



(30) 



Our model is a special case with F cxt (z,t) = —V'(z). 
Define P = mdz/dr and let (z,P) — (zo{t), Po(t)) be 
a solution of this equation with _Fld set to zero and 
(z, P) = (zo + Az, P + AP) be a linearized solution 
about (zq,Po) with Fld set to zero. Then Az and AP 
satisfy 

d 

di" 

" ' ' ' ' Az = B(t)Az. (32) 



: A: = mT 1 ^- zlf' 2 AP = A(t)AP , (31) 
dF cxt 



dt 



AP = 



dz 



Let (Az s (t), AP s (t)) be a solution of these equations 
satisfying (Az s (s), AP s (s)) = (0,1). Then (dz/dp) t = 
Azo(t). Hence Eq. l|2^ll becomes 



Szq 



dtF LD (t)Az (t). 



(33) 
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Now, the approximate solution (z + Sz,P + SP) of 
Eq. l(3T)|) with Fld 7^ obtained by perturbing (z ,P ) 
to first order in F LD satisfies 



d 

d_ 

dt 



Sz = A(t)5P, 
5P = B(t)Sz + F LD (t), 



(34) 
(35) 



where FhD{t) = Fld\z=z - The solution of these equa- 
tions with the initial condition lim t _ f _ 00 (<5z, 5P) = (0, 0) 
is given by 



Sz{t) = f dsF LB (s)Az s (t), 

J — OQ 

SP(t) = [ ds F LD (s)AP s (t) . 



Hence the classical position shift is 



Sz 



dtF w (t)Az t (0). 
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Thus, to establish that the position shift 5zq in QED at 
order e 2 in the limit h — > is equal to that in the Lorentz- 
Dirac theory, 5z, we only need to establish Az (t) = 
—Az t (0). In fact one can show in general that Az s (t) — 
—Az t (s) as follows. Note first that due to Eqs. 131fl and 
(|32|l the symplectic product defined by 

(Az« ApW\Az&, APW) = AzWapW-APWAzW 

(39) 

is time-independent for any pair of solutions 
(AzW, APW) and (Az< 2 \ AP^). Hence we have 

Az s (t)AP t (t) - AP s (t)Az t (t) 

= Az s (s)AP t (s) - AP s (s)Az t ( s ) . (40) 

Since Az s (s) = Az t (t) = and AP s (s) = AP t (t) = 1, 
we have Az s (t) = —Azt(s) as required. Hence 8zq = Sz. 

We have shown the equality of the classical and quan- 
tum position shifts at order e 2 in the limit h — ► for 
any system described by Eq. H3()|) provided that the pho- 
ton emission amplitude is given by Eq. flty) . We have 
derived Eq. I|18|l only for the case where F oxt (z,i) is t- 
independent although it appears likely to hold in a more 
general case. Let us verify that the classical position shift 
Sz with F e xt(z 7 t) — — V'(z) given by Eq. JSJ is indeed 
equal to Szq in Eq. Q29JI . though our general argument 
guarantees the agreement. To this end we need to cal- 
culate {dz/dp)t- Since the energy is conserved, we have 



J to 2 (dz/dr) 2 + m 2 + V(z) = y 'p 2 + m 2 . (41) 



Thus 



dz 
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By integrating this with the condition z = 
differentiating with respect to p, we find 



at t = and 



dz 
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where z\o — pj \Jp 2 



is the final velocity. The substi- 



tution of Eq. (|43|l in Eq. (|29|l gives the quantum position 
shift as follows: 



£|o 
m 



1 j dz 



(44) 

which is identical with the position shift Sz in the 
Lorentz-Dirac theory given by Eq. 

In this Letter we have reported that the h — * limit 
of the position shift due to radiation reaction in QED to 
lowest non-trivial order in e 2 is correctly reproduced by 
the Lorentz-Dirac theory for linear acceleration. It would 
be interesting to extend our result to three dimensional 
motion. Investigation of the one-loop correction to the 
potential is also important since this also affects the mo- 
tion of the particle. The corresponding correction to the 
position of the particle in a model with space-dependent 
mass term |6j was found to be of order ft -1 , overwhelming 
the effect from radiation reaction. However, in the more 
realistic model we have considered here, the one-loop cor- 
rection to the position appears to be at most of order h°. 
It would be very interesting to estimate this contribution 
and compare it with that from radiation reaction. 
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